Two stellar-mass black holes in the globular cluster M22 
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O 

Hundreds of stellar- mass black holes likely form in a typical globular star cluster, with all but 
one predicted to be ejected through dynamical interactions^. Some observational support 
for this idea is provided by the lack of X-ray-emitting binary stars comprising one black hole 
and one other star ("black-hole/X-ray binaries") in Milky Way globular clusters, even though 
many neutron-star/X-ray binaries are knowiP. Although a few black holes have been seen 
in globular clusters around other galaxies 5 6 , the masses of these cannot be determined, and 
some may be intermediate-mass black holes that form through exotic mechanisms 7 . Here 
we report the presence of two flat-spectrum radio sources in the Milky Way globular cluster 
M22, and we argue that these objects are black holes of stellar mass (each ~ 10-20 times more 
massive than the Sun) that are accreting matter. We find a high ratio of radio-to-X-ray flux 
for these black holes, consistent with the larger predicted masses of black holes in globular 
clusters compared to those outside 8 . The identification of two black holes in one cluster 

y—i shows that the ejection of black holes is not as efficient as predicted by most models^EEl, and 

we argue that M22 may contain a total population of ~ 5-100 black holes. The large core 

q radius of M22 could arise from heating produced by the black holes^. 

O 

y—i We have obtained very deep radio continuum images of the Milky Way globular cluster 

M22 (NGC 6656) with the Karl G. Jansky Very Large Array (VLA). The principal goal of the 
observations was to search for a possible central intermediate-mass black hole via synchrotron 

• £h emission from the accretion of intracluster gas; no central source was found 10 . However, we 

serendipitously detected two previously unknown radio continuum sources in the core of the cluster 
(Fig. 1). We term the sources M22-VLA1 and M22-VLA2. Both sources have flat radio spectra 
and are unresolved at our ~ 1" resolution. 



The core radius of M22 is uncommonly large for a Milky Way globular cluster: ~ 1.24 pcP^. 
These sources are well inside the cluster core, at projected radii of 0.4 pc and 0.25 pc for M22- 
VLA1 and M22-VLA2 respectively. The distance of the next source of comparable flux density is 
far outside the core, at a projected radius of 2.4 pc. 
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These sources have no counterparts in shallow archival Chandra X-ray imaging. Based on 
these non-detections, the sources are constrained to have Lx ~ 2.2 x 10 30 erg s _1 over 3-9 keV 
at the distance of M22r^. The radio luminosities of the sources are Lr ~ 6 x 10 27 erg s _1 at 8.4 
GHz, assuming flat spectra. Therefore, if the sources are not variable, the limit of radio to X-ray 
luminosity is: log L R /L X ~ — 2.6. 

The radio luminosity, Lr/Lx ratio, and central location of the sources place significant 
constraints on their nature. The most likely explanation is that both sources are accreting stellar- 
mass black holes in M22. Other possibilities, all of which we consider unlikely, are discussed in 
the Supplementary Information. These objects are the first strong candidates for stellar-mass black 
holes in any Milky Way globular cluster, and the first stellar-mass black holes to be discovered 
through radio emission rather than via X-rays.^ 

The radio emission implies that the black holes are actively accreting, and the flat radio spec- 
tra are consistent with relatively low accretion rates^ ( <, 2-3% of the Eddington rate). Because 
globular clusters have modest amounts of interstellar gas, it is very unlikely that the radio luminos- 
ity can be explained by Bondi accretion. Thus the objects cannot be black-hole/black-hole binaries, 
and instead are probably in binary systems with Roche lobe-overflowing companions. Stellar-mass 
black holes (~ 5-100 M ) offer the best explanation for the presence of multiple sources close to 
the cluster center; objects more massive than the average cluster star will sink to the center because 
of mass segregation. 

To look for optical counterparts to the radio sources, we used archival Hubble Space Tele- 
scope (HST) imaging of M22, for which photometric catalogs are available 15 . Fig. 2 shows that 
M22-VLA1 is a close match (0.05") to a moderately low-mass (~ 0.34 M ) main sequence M 
dwarf in M22 as inferred from standard stellar isochrones (see Supplementary Information for 
more details). M22-VLA2 is 0.17" from the nearest detected star, which is a ~ 0.62 M main 
sequence star. Considering the distribution of stars in the inner 30" of the cluster, the probability 
of a chance coincidence as close as for M22-VLA1 is only 2%; for M22-VLA2 it is 26%. Thus we 
consider the optical association for source M22-VLA1 suggestive, but that for M22-VLA2 uncer- 
tain. However, for the case of M22-VLA1, there is an additional complication: since the average 
stellar mass in the core is greater than that of the putative companion, the low-mass main sequence 
star would likely be exchanged out of the binary in a three-body interaction with another staP. On 
the other hand, because of the low central density of M22 1 1 (< 10 4 M Q pc~ 3 ), a binary with a low- 
mass companion may survive longer than in a typical globular cluster. Nonetheless, it is possible 
that both radio sources are associated with low luminosity objects below the detection limit of the 
HST data, such as white dwarfs. 

Stellar-mass black holes with accretion rates below ~ 2% of the Eddington rate 16 (in the 
so-called low/hard state) follow an empirical correlation between radio and X-ray luminosity 
with a scatter of ~ a factor of twcP-^. Fig. 3 shows this correlation with the M22 data overplot- 
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ted. The radio-X-ray relation predicts an X-ray luminosity of 10 31 -10 32 erg sec -1 for this radio 
luminosity-^EI], above the completeness limit of the archival Chandra data. There are several plau- 
sible explanations for this discrepancy. First, there is the possibility of variability. The X-ray data 
were taken in in 2005, six years earlier than the radio data. Field stellar-mass black holes in the 
low/hard state show substantial (typically a factor of 2-10) variability in both radio and X-rays^lMJ 
Therefore, concurrent radio and X-ray data are necessary for precise constraints on L R /L X - We 
found marginal evidence for radio variability in M22-VLA2 on the timescale of a week; more de- 
tails can be found in the Supplementary Information. Another plausible explanation is that there is 
larger scatter in the radio-X-ray correlation at very low accretion rates. Only a single known black 
hole binary has a measured radio luminosity as faint as our sources^, and there is evidence that 
some stellar-mass black holes with low X-ray luminosities may not fall on the correlation^. 

An intriguing possibility is that these sources have high values of L R /L X because they are 
more massive than typical stellar-mass black holes in the field. The radio-X-ray correlation for 
stellar-mass black holes is a special case of a "fundamental plane" for black hole accretion in 
the low/hard state that includes the black hole mass as a third parameter 22 . In this relation, more 
massive black holes have larger values of Lr/Lx- If our sources have masses of ~ 15-20 M Q 
rather than the 5-10 M Q typical of field stellar-mass black holes 23 , then their X-ray luminosities 
should be lower than predicted by the correlation in Fig. 3 by a factor of ~ 2-3. It is reasonable 
to expect that black holes in globular clusters will be more massive than those in the field. Field 
black holes with measured dynamical masses are all in binary systems and were probably affected 
by mass transfer during a common envelope stage that reduced the mass of the resulting black 
holes 24 . This need not be the case in globular clusters, since black holes can form as single objects 
or in wide binaries, and then be exchanged into pre-existing binaries or tidally capture companions 
due to the high stellar densities^. Globular cluster black holes also form at lower metallicity than 
in the field, leading to less mass loss from the progenitor and thus more massive remnants 8 . 

As mentioned above, the location of stars in a cluster also gives information about their 
masses. Stellar-mass black holes will mass-segregate to the core of the cluster. This process can be 
used to roughly estimate their masses by assuming thermalization, for which this relation holds 1 : 
mBH/ m ± = ( r c/ r BH) 2 , where m BH and r BH are the characteristic black hole mass and radius, 
is the typical stellar mass, and r c is the core radius. Assuming = 1 M in the segregated cluster 
core and taking the observed values of r c = 1.24 pc and r BH = 0.33 pc, we estimate m BH ~ 15 



The existence of black holes in a low-density globular cluster such as M22 constrains the 
magnitude of initial velocity kicks received by the black holes at birth. The current central escape 
velocity of M22 11 is ~ 34 km s _1 . This value may have been higher in the past, due to a larger 
cluster mass and a more compact structure. Nonetheless, the retention of two black holes in a 
globular cluster with a modest escape velocity implies that the black holes could not have received 
large natal kicks. Large kicks are inferred for some stellar-mass black holes in the field 2 ^. Low kick 
velocities can originate from supernovae if the black hole mass is large, or if the black holes form 
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from direct collapse with no supernovae. In either case, higher black holes masses are favored. 

The presence of black holes in a globular cluster can lead to an expansion of the core radius 
through interactions between black holes and stars. This could explain why M22 has the fifth- 
largest core radius among luminous ( £ 2 x 1O 5 L ) Milky Way globular clusters^. Additional 
discussion can be found in the Supplementary Information. 

Most theoretical models in the literature predict that only a single black hole (or black- 
hole/black-hole binary) will survive the dynamical processes by which black holes mass-segregate 
to the cluster center, form an unstable subcluster, and evaporateU^I. In some cases more than one 
black hole may temporarily survive for an additional black hole relaxation time (< 1 Gyr), if the 
extra black holes are kicked into orbits outside the coreP. Additional discussion can be found in 
the Supplementary Information. 

In contrast to these theoretical predictions, M22 contains more than one black hole. In fact, 
it is possible that more than two black holes are present in M22, either as single black holes or in 
binary systems that are not undergoing observable mass transfer. Under the uncertain assumption 
that both of the M22 sources are black hole-white dwarf binaries, published calculations can be 
used to estimate the fraction of surviving black holes that are actively accreting in present-day 
globular clusters^. Over 10 Gyr, 2-A0% of black holes are expected to become members of binary 
systems with observable accretion. Our two observed sources thus suggest a total population of 
~ 5-100 black holes in M22. 
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Figure 1: VLA radio continuum image of the core of the globular cluster M22. The two bright 
circled objects are the sources identified as stellar-mass black holes, M22-VLA1 and M22-VLA2. 
These sources have flux densities of 55-58 /iJy at 5.9 GHz. We obtained the data in two separate 1 
GHz basebands centered at 5 and 6.75 GHz, allowing a measurement of the spectral index of the 
radio emission between these frequencies. Both sources have flat radio spectra, with a = — 0.2, 
assuming S u = v a . The faint circled object is a known millisecond pulsar^. A red cross marks 
the photometric cluster center. 20" corresponds to approximately 0.3 pc at the distance of M22. 
The apparent elongation of the sources is due entirely to the elongated synthesized beam; all three 
circled sources are unresolved. North is up and east is to the left in this image. 
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Figure 2: Optical images of M22 and the candidate companion stars to the radio sources, (a) 
Ground-based image that shows the approximate location of the sources in the context of the 
star cluster, (b) and (c) show the zoomed-in location of the radio sources on an archival high- 
resolution Hubble Space Telescope/ 'Advanced Camera for Surveys F814W image. These circles 
have radii of 0.3" for clarity; the uncertainty in the astrometric matching of the optical and radio 
data is < 0.1". The image orientation is as in Fig. 1. (Image credit for (a): Doug Matthews/ Adam 
Block/NOAO/AURA/NSF.) 
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Figure 3: Radio-X-ray correlation for stellar-mass black holes. The M22 sources have properties 
more consistent with black holes than neutron stars or white dwarfs. Filled squares represent si- 
multaneous radio and X-ray data; unfilled squares are non- simultaneous measurements for which 
variability may affect their positions. Upper limits are also shown. Some objects have multiple 
measurements plotted that represent different phases of accretion. The open red circle represents 
both M22-VLA1 and M22-VLA2, which have very similar luminosities. The dotted black line 
represents the published correlationP' L R oc L^ 58 , normalized by a least- squares fit to the simul- 
taneous detections with L x < 2 x 10 34 erg s _1 . The dashed and dot-dashed blue lines show two 
possible radio-X-ray correlations for accreting neutron stars; this relation is poorly constrained by 
observations 28 . The solid green line shows the maximum radio continuum luminosity observed for 
accreting white dwarfs 29 . Neither neutron stars nor white dwarfs have properties consistent with 
the M22 radio sources. More information can be found in the Supplementary Information. 
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Supplementary Information 



1 Data Acquisition and Reduction 

We observed M22 with the VLA in May 201 1 as part of the program 10C- 109 (P.I. L. Chomiuk). 
Ten hours were spent observing the cluster, split amongst four 2.5 hour blocks, for a total of 7.5 
hours on each source. We observed in C band with 2 GHz total bandwidth and four polarization 
products. Of the two basebands of 1 GHz width, one was centered at 5.0 GHz and the other at 
6.75 GHz. The array was in BnA configuration, which is the normal B configuration with an 
extended northern arm to permit improved observations of southerly sources. This configuration 
gave a resolution of 1.53" x 0.81" at our mean frequency of 5.9 GHz. The field of view (full- 
width at half power) of the VLA at this frequency is ~ 7.7' in diameter, compared to the 
half-mass diameter of 6.1' for M22P. 

We observed J1820-2528 as the secondary phase calibrator and J1407+2827 as the polar- 
ization leakage calibrator. 3C286 was used as an absolute flux density, bandpass, and polar- 
ization angle calibrator. The data were reduced using standard routines in AIPS. Weights were 
applied using TYAPL, and each 2.5-hour block was edited for bad data and interference, fol- 
lowed by calibration. For each individual calibrated baseband, the data were concatenated in 
the uv plane and then imaged. Figure 1 is a deep coadded image of both basebands, obtained 
by smoothing the 6.75 GHz baseband to the resolution of the 5.0 GHz basebands and averaging 
these together. The rms sensitivity of this coadded image is 1.5 /Jy beam -1 , one of the deepest 
radio continuum images ever published. 

We measure flux densities on the individual image for each baseband. The measured flux 
densities for the two sources are: M22-VLA1: 56.3 ± 3.9 /Jy (5.0 GHz), 60.1 ± 3.6 /Jy (6.75 
GHz); M22-VLA2: 54.6 ± 4.3 /Jy (5.0 GHz), 54.5 ± 3.5 /Jy (6.75 GHz). If we assume the 
flux densities follow a power law of the form S u = v a , a = 0.2±0.3 and 0.0 ± 0.3 for sources 
M22-VLA1 and M22-VLA2 respectively. We set upper limits of ~ 10% on both the linear and 
circular polarization of the sources at C band. 

Matching of the VLA and HST astrometric systems was required to search for opti- 
cal counterparts to the radio sources. We fixed the astrometric zeropoint of the F814W (J- 
equivalent) HST image to the VLA data through a match between an obvious background 
galaxy in the HST image and a resolved continuum source in the VLA images. This map- 
ping was confirmed to a precision < 0.1" with the X-ray image, which also has sources in 
common with each of the HST and VLA data. The J2000 positions of the sources are, in the 
reference frame of the VLA images: M22-VLA1 (18:36:25.825 -23:54:13.75), M22-VLA2 
(18:36:23.824 -23:54:33.50). The formal uncertainties in these positions are ~ 0.02 - 0.03". 

2 Calculation of X-ray Upper Limits 

The archival Chandra ACIS-S imaging has a total integration time of 16 kiloseconds. If we 
assume a completeness limit of 8 counts, then at the distance of M22 (3.2 ± 0.3 kpc 12 ) these 
data are sensitive to a luminosity of L x ~ 2.2 x 10 30 erg s _1 over the range 3-9 keV assuming 



a foreground N H = 1.7 x 10 21 cm 2 and a photon index of 2. If instead taken over the Chandra 
energy band of 0.5-8 keV, the limit is ~ 5.6 x 10 30 erg s _1 . 

3 Photometry of Candidate Optical Counterparts 

The optical photometry of the candidate counterparts 15 , corrected for differential reddening^, 
is: FQ0QW = 21.77±0.11; (FQOQW -F8UW) = 1.04±0.11 (M22-VLA1) and FQOQWq = 
18.62 ± 0.01; (FQOQW - F8UW) = 0.59 ± 0.02 (M22-VLA2). The listed photometric 
uncertainties do not include the uncertainty in the extinction. For a distance of 3.2 kpc, the 
inferred absolute magnitudes are M F8UW = 8.21 and 5.51 respectively. These values were 
converted into main sequence masses using a standard stellar isochrone for these parameters^ 
13 Gyr, [Fe/H] = -1.7, [a/Fe] = +0.2. 

Compared to the ridge line of main sequence stars in M22, the candidate counterpart for 
M22-VLA1 lies 0.08 ± 0.11 mag blueward in (FQOQW - F8UW). A blue excess could be 
taken as evidence for hot emission from a putative accretion disk, but given the large error in 
the color the difference is not statistically significant. 

4 Other Candidates for the Sources 

We have discussed the classification of the radio continuum sources as stellar-mass black holes 
in the main text. We believe an unlikely (but possible) alternative classification is as background 
sources, which we consider first. Then we discuss other possible classifications that we believe 
are very unlikely or can be ruled out. Future very long baseline interferometry can be used to 
measure the proper motions of M22 VLA-1 and VLA-2 and definitively confirm or refute their 
association with M22. 

Background sources: Very little information exists about the radio properties of background 
objects below 100 /iJy in flux density, especially when spectral index information is included. 
Using the limited deep background source counts from the literature^, we find that ~ 0.2-0.4 
background objects with flat spectra and flux densities of 30-100 /iJy would be expected in the 
1.3' radius core of M22. In the central 30", < 0.1 source would be expected. These calculations 
do not account for the lack of optical counterpart, which reduces the expected counts by a factor 
^2 — 3 (see below). For these reasons, we believe that it is very unlikely, but not impossible, 
that M22-VLA1 and M22-VLA2 are both background sources. 

Background galaxies can have flat radio spectra through two scenarios: active galactic nu- 
clei with partially self-absorbed jets, or star-forming galaxies in which the synchrotron emission 
is "diluted" with thermal emission to produce a flat composite spectrum. Active galactic nuclei 
are expected to dominate at the flux densities of our sources. Most background radio continuum 
sources in the observed range of flux densities are found to have optical galaxy counterparts 
brighter than the HST completeness limit of / ~ This is especially true of sources 

with flat radio spectra, many of which are associated with relatively bright early-type galaxies 
at z < f2. There are no candidate optical background galaxies that are plausibly associated 
with either source. By contrast, in the outer regions of our M22 VLA image there is a canonical 
background galaxy: it has a flux density of ~ 100/iJy at 5.9 GHz and is resolved in the VLA 



data, it has a spectral slope a = —0.8, and has an obvious distant early-type galaxy as an optical 
counterpart in the HST images. 

Two published studies with multi-band radio continuum imaging can be used to assess the 
occurrence of background dopplegangers to the M22 sources. The first uses deep VLA 1.4 and 
8.4 GHz data and HST imaging^. Of a catalog of 65 sources with optical data, only two have 
inferred 5.9 GHz flux densities in the range 30 to 100 /xJy and a = —0.2 to +0.2, and both 
have bright (/ < 20) early-type galaxy counterparts. Another study used VLA 1.4 and 4.8 GHz 
data in the Chandra Deep Field South 32 . Here there are 12 sources that satisfy these flux density 
and slope ranges, though some of the spectral slopes have large uncertainties. Eight of 12 are 
bright enough that they would have been detectable in the HST data; two others are borderline. 
From these (limited) comparisons we conclude that the lack of background optical counterparts 
reduces the expected rate of radio sources with the observed properties by at least a factor of 
2-3. 

Pulsars/Supernova remnants: Millisecond pulsars are commonly found in globular clusters. 
However, these have universally steep radio spectra (a £ — 1), inconsistent with the flat 
measured spectra. In rare cases pulsars drive a wind that interacts with ambient material, termed 
a pulsar wind nebula. These objects have flat measured spectra and Lr/Lx consistent with the 
constraints for our sources. However, they also generally have L x > 10 34 erg/s, a high degree 
of radio polarization, large sizes, are short lived, and are found in regions with dense gas, 
all inconsistent with observations 38 . Supernova remnants have similar properties and are also 
excluded. 

Accreting compact objects: A number of candidates involve accreting neutron stars or white 
dwarfs. For most companions these systems will have much lower Lr/Lx or Lr than accreting 
black holes, so are not consistent with the observed constraints on our sourcesS^MEQ Figure 3 
shows this graphically. The bl ack hole data in this figure are compiled from recent papers on 
stellar- mass black holed ^l 42 !! 3 ]. 

A contrived exception to this argument is if the sources are very strongly variable; the radio 
luminosities could be consistent with flares from accreting neutron stars, which then might have 
Lx just below that which would have set off all-sky X-ray monitoring satellites. However, this 
scenario would require the coincidence of both neutron stars flaring simultaneously. In addition, 
they would need to be abnormally faint in quiescence^, since the typical L x for quiescent 
accreting neutron stars is > 10 32 erg s _1 . 

Symbiotic stars, which are white dwarfs accreting from luminous red giants, are ruled out 
by the optical source matching. 

Planetary Nebulae: Planetary nebulae can emit optically-thin thermal emission at radio wave- 
lengths, and would therefore show a flat spectral index consistent with our sources. However, 
planetary nebulae are also bright [O III] 5007A emission line sources. No [O III] nebulosity is 
observed around either source in archival HST/WFPC2 images obtained in the narrow F502N 
filter. By contrast, the known M22 planetary nebula IRAS 18333-2357 is clearly detected in this 
HST [O III] image, and is not detected in our VLA radio continuum image. 



Foreground ultracool dwarfs: The last remaining possibility is that the sources are nearby 
(< 100 pc), cool (late M to early L) active dwarf stars in the foreground, which satisfy the 
Lr/Lx constraints and which can have flat radio spectra 45 . The stars cannot be too nearby 
( ~ 50 pc) or else they would be detectable in the HST data, so they would need to be in a 
relatively narrow range of ~ 50 — 100 pc, the larger end of which would require strong flares 
(and be inconsistent with the lack of strong radio variability of our sources). Some ultracool 
dwarfs also show strong circular polarization 46 , which we do not observe. Cool dwarfs will 
be more luminous in the infrared, so we have consulted archival K data from the Vista VVV 
survey 47 . These data rule out a cool dwarf identification with source M22-VLA1. The constraint 
for source M22-VLA2 is less strong as it is located near the diffraction spike from a bright star. 
Thus we cannot fully discount this possibility for M22-VLA2, but consider it to be unlikely 
because of the rarity of these sources. 

5 Variability 

We searched for variability in the radio sources by imaging each 2.5-hr data block separately 
and measuring their flux densities at each epoch. Fig. SI shows the results. Over the week of 
observations, M22-VLA1 appears constant within the uncertainties, while the flux density of 
M22-VLA2 appears to decrease monotonically by about 20%. This overall trend is significant 
at the ~ 2.6a level. Thus, the evidence for variability is only suggestive with the present data, 
and additional observations are necessary to determine if significant variability is present. 

6 The Large Core Radius of M22 

The core radius of M22 is 1.24tgo| pc, which is fifth-largest among luminous (^2x 1O 5 L ) 
Milky Way globular clusters (it is third-largest, next to M14 and M53, if the anomalous clusters 
co Cen and NGC 2419 are excluded) 11 . A number of processes can lead to an increase in the 
core radius for globular clusters that contain black holes, and could help explain the large core 
radius of M22. The process dominant in some theoretical simulations 9 is the formation of one or 
more black hole-black hole binaries in the cluster core. Three-body interactions between one of 
these binaries and another black hole can eject the black hole from the core (but not the cluster 
entirely), and as this black hole sinks back to the center through dynamical friction it heats stars, 
removing them from the core. It is not clear what the necessary minimum population of black 
holes is for this process to be effective^'. 

It is also possible that binaries containing a black hole and another object, such as a main 
sequence star or white dwarf, can also heat the core through three-body scattering interactions. 
Black hole binaries with Roche lobe-overflowing counterparts can have binding energies (~ 
10 48 -10 50 erg s _1 ) comparable to the total binding energy of all of the single stars in the core 
of a globular cluster. However, these interactions may instead lead to stellar mergers^. 

A possibility frequently discussed in the literature is that stellar interactions with an 
intermediate-mass black hole (> 100 M ) can also lead to heating and core expansion 50 51 . 
Deep radio data for M22^have ruled out central intermediate-mass black holes with masses 
> 360 M Q (3a) but less massive black holes could still be present. These would have a mass 



ratio < 6 x 10 4 compared to the cluster mass and so it is unclear whether they could be solely 
responsible for a significant increase in the core radius, but some contribution is possible. 

Finally, it is possible that the large core radius of M22 is unconnected to the presence of 
black holes. It has been proposed that the dynamical evolution of most globular clusters is less 
advanced than generally supposed, and that many clusters are still in an initial, slow phase of 
core contractior l 52 ! 53 !. In this case, the unusually large core radius of M22 might just reflect the 
initial conditions of the formation of the cluster. 

7 The Origin of Multiple Black Holes in M22 

Here we mention a few possible explanations for the presence of multiple black holes in M22. 
The discussion in the main text emphasized theoretical work suggesting that only one black hole 
or black hole-black hole binary will generally remain through long-term dynamical evolution 
of a globular cluster. A contrasting view is drawn from simulations in which the presence of 
black holes causes significant core expansion 9 . This expansion leads to a large reduction in the 
stellar density in the core and slows the black hole interaction rate. Under these circumstances, 
a significant fraction of the initial complement of black holes remain bound over timescales 
of several Gyr. These simulations were designed to match Magellanic Cloud star clusters with 
lower masses and larger core radii than M22 and so have central densities much lower (by a 
factor > 10 3 ) than M22 at the present day. Comparable simulations tailored to the properties of 
M22 can help constrain the formation of M22-VLA1 and M22-VLA2. 

A speculative hypothesis for the existence of multiple black holes arises from the observa- 
tion that M22 has modest spread in [Fe/H] 54 . [Fe/H] dispersions are common among massive 
globular clusters and may be due to self-enrichmenP^ it has also been proposed (for M22 and 
other massive clusters) that an [Fe/H] spread is due to the merger of two or more globular 
clusters that were each homogeneous in [Fe/H], but with different mean levels of enrichmenP^. 
Such mergers could occur either in dwarf galaxies with low velocity dispersions or if the clusters 
were in binary systems. If these putative mergers happened after the formation and dynamical 
evolution of black hole subsystems in globular clusters, then, even if each cluster originally had 
only one remaining black hole, multiple black holes could be present after the merger. 
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Figure SI: Time series of flux densities for M22 radio sources. M22-VLA1 (open circles) is 
constant over the week of observations, while M22-VLA2 (filled circles) decreases over the same 
period at a significance of 2.6(7. The error bars are standard deviations of the measured flux densi- 
ties. 



